Radical-free dynamic nuclear polarization using electronic defects in silicon 
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Direct dynamic nuclear polarization of X H nuclei in frozen water and water-ethanol mixtures is 
demonstrated using silicon nanoparticles as the polarizing agent. Electron spins at dangling-bond 
sites near the silicon surface are identified as the source of the nuclear hyperpolarization. This novel 
polarization method open new avenues for the fabrication of surface engineered nanostructures to 
create high nuclear-spin polarized solutions without introducing contaminating radicals, and for the 
study of molecules adsorbed onto surfaces. 



Dynamic nuclear polarization (DNP) has recently 
emerged as a powerful technique for improving the sensi- 
tivity of nuclear magnetic resonance (NMR) experiments, 
leading to new insights in materials characterization p]- 
4J, complex biomaterials [5-8J, and diagnostic medicine 
using magnetic resonance imaging (MRI) [9-12J. By ap- 
plying microwave (/iW) irradiation at or near the electron 
Larmor frequency, a sizable nuclear spin polarization can 
be generated via the transfer of electron spin polariza- 
tion to nearby nuclear spins. The effect can be enhanced 
by operating at low temperatures or high magnetic fields, 
where the equilibrium electron polarization can approach 
unity. The electrons used in the polarization process may 
originate from free radicals or paramagnetic metal ions 
introduced into the system, or be intrinsic to the material 
being polarized. Both the rate and saturation magnitude 
of nuclear polarization depend strongly on the concentra- 
tion of unpaired electrons. This is because spin transport 
via nuclear spin diffusion is slow [13j [14] and competing 
relaxation processes means that a high radical concentra- 
tion is required to achieve large polarizations on experi- 
mentally relevant time scales. 

Most biological applications of DNP involve systems 
that have no intrinsic unpaired electrons for polariza- 
tion, and so a significant effort has been undertaken 
to develop a range of free radicals fTSHlT] that can be 
dissolved in a solvent (such as water, ethanol, or glyc- 
erol) along with the biomolecular substrate to be polar- 
ized before freezing for DNP. However, the presence of 
these radicals in the polarized substrate after polariza- 
tion reduces the spin lattice relaxation time (Ti) of the 
polarized substrate p~8| [19], interferes with the chemi- 
cal environment of molecules under study [20J and, for 
biological experiments, poses significant toxicity issues. 
To reduce contamination of the sample by residual radi- 
cals, filtering [21J and radical scavenging [18J undertaken 
post-disolution have been attempted with limited suc- 
cess. The radicals may be bound into a porous organic 



matrix [22, 23J, though filtration is still necessary. Al- 
ternatively, the source of polarization may be physically 
separated from the material to be polarized. It has been 
shown that the large nuclear spin polarizations present in 
optically pumped noble gases and semiconductors can be 
transferred to a variety of materials using a Hartmann- 
Hahn cross-polarization sequence [24-27J or, in the case 
of gases, via a direct nuclear dipole-dipole interaction 
[28] . Surface defects in porous carbon chars have also 
been shown to polarize 1 H nuclei in solution at room 
temperature via a contact hyperfine interaction |29j . 

In this Letter, we demonstrate that naturally occurring 
paramagnetic defects at the surface of silicon nanopar- 
ticles (SiNP) can directly polarize 1 H nuclei in a bulk 
frozen matrix near the surface of the particle as well as 
the 29 Si nuclei in the particles. Two distinct nuclear spin 
baths are observed, one consisting of nuclei polarized di- 
rectly by surface electrons, and a second consisting of 
nuclei polarized indirectly by nuclear spin diffusion. By 
comparing natural abundance and isotopically enriched 
samples, we find that 1 H polarization is strongly affected 
by the isotopic concentration of 29 Si nuclei inside the 
particles. This technique allows uncontaminated hyper- 
polarized solutions to be generated at room temperature, 
benefitting applications from polarized targets to molec- 
ular imaging. 

Samples consisted of natural abundance ( 29 Si = 4.8%) 
polycrystalline powder (Alfa Aesar, d = 3 /im), and iso- 
topically enriched ( 29 Si = 91.4%) [30] silicon particles 
of diameter d — 200 nm made by ball milling bulk sili- 
con and then size-separated by centrifugal sedimentation 
[31] , Particle sizes were confirmed by scanning electron 
microscopy (SEM). Suspensions of these particles in wa- 
ter (1:1 Si:H20 by weight), and ethanol and water (10:9:1 
Si:EtOH:H20 by weight) were prepared by sonication and 
then degassed with a freeze-pump-thaw cycle before be- 
ing loaded into the cryostat for measurement. 

DNP experiments were performed at 4.2 K using a con- 
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Figure 1. Microwave frequency dependence of the (a) 
X H NMR signal (closed circles) and (b) 29 Si NMR signal 
from a sample of natural isotopic abundance SiNP suspended 
in frozen water, measured at 4.2 K. No change in the 1 H NMR 
signal was observed for the same sample without the SiNP 
(open circles), and no 1 H signal was observed from a dry sam- 
ple of SiNP. Upper insets in (a) and (b): X H and 29 Si spectra 
for /xW irradiation frequency of 65.875 GHz. Lower inset 
in (b): Room temperature ESR signal of natural abundance 
SiNP. 



tinuous flow cryostat in an applied field, B, of 2.35 T 
(cj 29Si = 20 MHz, c^ H = 100 MHz). The sample was 
housed in a glass capillary tube in direct contact with 
the flowing helium vapor. A Bruker NMR spectrometer 
with a home-built probe and solenoidal coil [32] were used 
for NMR detection. Microwave irradiation was provided 
through a 90 mW Gunn diode source (Millitech), and the 
/iWs were coupled from room temperature to the sam- 
ple via a mm waveguide and horn antenna. The nuclear 
polarization was measured with a saturation recovery se- 
quence ( 29 Si) (7r/2) 16 — t — (tt/2) x or saturation recovery 
solid echo sequence ( 1 H) (tt/2) —t — (7r/2) x —td — (7r/2) Y 
due to the rapid decay of the X H signal, corresponding to 
a ~ 100 kHz X H dipolar linewidth. 

Figure 1(a) shows the /iW frequency dependence of the 
X H NMR signal for solid H2O samples, with and with- 
out suspended SiNP (natural isotopic abundance), for 
a fiW irradiation time of 120 s. For the sample with 
SiNP, the X H nuclear polarization depended strongly on 
the applied /iW frequency, while no change in X H polar- 
ization was seen for the sample without SiNP. A similar 



Figure 2. (a) Schematic model showing the location of DNP- 
active electron spins at and near the silicon-silicon oxide in- 
terface. Arrows indicate the directions of nuclear polarization 
flow away from the surface of the particle, (b) Thermody- 
namic model (see text) for the polarization process of both 
the *H and 29 Si nuclear spin baths. 29 Si and 1 H nuclei near 
the defect electrons are polarized directly through a dipolar 
interaction on time scales rsn and thi- Nuclei far away from 
the surface are then polarized by nuclear spin diffusion on 
time scales rs%2 and thi that are dependent on the spatial 
separation of the nuclear spins. 



dependence on jaW frequency of the 29 Si NMR signal for 
the same sample of SiNP suspended in frozen water (Fig. 
1(b)) after 600 s of /iW irradiation. No X H signal was ob- 
served from a sample of dry SiNP. Representative X H and 
29 Si spectra at a fiW frequency of 65.875 GHz are shown 
in the insets [33J. We note that no change in the room 
temperature X H spectrum was seen in a similarly pre- 
pared sample after repeated measurements over a period 
of a week. This suggests that chemical reactions at the 
silicon/water interface leading to Si-H bond formation is 
insignificant, and that the 1 H NMR signal results from 
X H nuclei outside of the particle rather than absorbed 
into the oxide or particle itself. 

Both X H and 29 Si polarizations show inversions at 
the same /iW frequency, ujo = 65.95 GHz, indicating 
that polarization of both nuclear species originates from 
electrons with the same g- factor. Room temperature 
X-band electron spin resonance (ESR) (lower inset of 
Fig. 1 (b)) of the particles shows a single broad peak 
at B = 335.9 mT and frequency 9.444 GHz, correspond- 
ing to g = 2.006. This g-factor is consistent with the 
orientation-averaged signal from unpaired electrons at P5 
defects, which are known to occur at the silicon-silicon 
dioxide interface |34] . To confirm the role of these defects 
in the polarization process we prepared a similar sample 
where the surface oxide was removed by a hydrofluoric 
acid etch in an inert atmosphere and suspended in glyc- 
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Figure 3. Time evolution of the 1 H nuclear polarization with 
(solid shapes) and without (open shapes) applied /xW irra- 
diation for frozen solutions of SiNP : H 2 (4.7% 29 Si) and 
SiNP : EtOH : H 2 (4.7% 29 Si, 91.4% 29 Si). The enhance- 
ments are scaled relative to the 4.2K equilibrium polariza- 
tion. Fits are exponential (dashed lines) or biexponential 
(solid lines), see text. Error bars are shown when greater 
than the data symbol. 

erol. This sample showed no ESR signal and no 1 H or 
29 Si DNP enhancement under similar conditions. Addi- 
tionally, the sign of polarization enhancement is the same 
for both nuclear species for frequencies either side of cjo, 
despite the gyromagnetic ratios of the X H and 29 Si nu- 
clei being opposite in sign. From this we can determine 
that the 1 H nuclei are being polarized through a direct 
interaction with the electron dipolar spin bath, rather 
than via a dipolar interaction with polarized 29 Si nuclei 
close to the surface. The direct interaction scales oc 7 2 
[37J, independent of the sign of 7 n , the nuclear gyromag- 
netic ratio. If instead the 1 H polarization was mediated 
by a dipolar interaction with the polarized 29 Si nuclei, 
opposite nuclear polarizations would be generated in the 
two nuclear species for the same polarizing frequency. 

A schematic model of the system is shown in Fig. 2. 
The anisotropic g factor of the P5 defect (g = 2.002 — 
2.010) [34J, together with the high density of unpaired 
electrons at the sample surface (~ 10 13 cm -2 ) [35J yields a 
broad spectrum of electron-electron dipolar interactions 
with energies at or near the Lamor frequency of each 
nuclear species. This results in a broad fiW frequency 
response of both X H and 29 Si nuclei (Fig. 1), consistent 
with previous DNP studies of silicon microparticles |3J. 

The time evolution of the X H polarization of vari- 
ous frozen suspensions of natural abundance and 29 Si- 
enriched SiNP in frozen water and water-ethanol mix- 
tures at /iW frequency of 65.875 GHz is shown in 
Fig. 3. All data are normalized to the 1 H equilib- 
rium polarization (P#) at 4.2 K without /iW irradia- 
tion. The enhancement curves with /iW irradiation dis- 
play biexponential behavior, well fit by the form P = 
Poo (1 - (Ae-^^ 1 + Pe-^ 0l / T2 )) (solid lines). This is 
in contrast to the time evolution of the X H NMR sig- 



Sample Ti(s) r H i(s) r H 2(s) 

Si:H 2 (4.7% 29 Si) 1113 ± 96 120 ± 14 2304 ± 307 
Si:EtOH:H 2 (4.7% 29 Si) 4.3 ± 0.4 4.6 ± 0.8 217 ± 58 
Si:EtOH:H 2 (91.4% 29 Si)| | 16 ± 3 | 97 ± 17 



Table I. 1 H polarization times fitting to experimental data in 
Fig. 3. T\ is the characteristic time taken for the 1 H nuclei 
to return to the equilibrium polarization without /iW irradia- 
tion, thi and TH2 are fast and slow polarization times which 
form two additive components of the polarization process cor- 
responding to polarization dominated by direct DNP and spin 
diffusion respectively. 

nals without /iW irradiation, which are better described 
by a single exponential, P = Pg (l — e _£pol / Tl ) (dashed 
lines). Best fit values for these characteristic times are 
given in Table 1. 

Note in Fig. 3 that the long-time induced polarization 
in the water-ethanol sample is suppressed relative to the 
water sample. We attribute this to fast nuclear relax- 
ation caused by the rotating methyl group in ethanol, 
which creates an oscillating magnetic field with a charac- 
teristic time constant of a few seconds even at cryogenic 
temperatures (methyl rotamer effect |36j). Increasing the 
29 Si isotopic concentration from 4.7% to 91.4% and re- 
ducing the particle size from 3 /im to 200 nm reduces 
the polarization rates and the total 1 H polarization en- 
hancement by a factor of two, from 1.5 to 1.25 times the 
Boltzmann polarization at 4.2 K. This is particularly no- 
table given the higher number of polarization sites arising 
from surface area to volume ratio that is 15 times larger 
in these samples. Despite this significant increase in the 
availability of directly polarization sites, the higher con- 
centration of 29 Si nuclear spins in the NP are able to 
absorb a greater total angular momentum transfer from 
the surface electrons compared to the natural abundance 
samples. This results in a slower direct polarization rate 
of the X H nuclear spins close to the surface. However, the 
relaxation due to the methyl rotamer effect and other 
electronic defects within the frozen matrix remains the 
same as for the natural abundance samples, and so com- 
petes more strongly with polarization by spin diffusion 
and hence causes the 1 H polarization to saturate earlier 
and at a lower absolute value. 

It is possible to estimate the size of the region of direct 
polarization for each nuclear species. Upon /iW irradia- 
tion of the electron spin bath, the polarization of nuclei 
within a radius r < j3 is through a direct dipolar interac- 
tion with a pair of electrons at a rate r^ 1 = T/r 6 , where 
T = KCG(ujq — oj*), K = 7T7 e P 2 r/2 is the saturation 
strength of the microwave field, 

10 B 2 1 + B 2 7 2T 2 V { } 

is a constant describing the nature of the electron- nuclear 
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dipolar interaction, and G(u — cj*) is the electron line- 
shape function [37J. Here B\ is the strength of the 
fiW magnetic field, r is the correlation time of the elec- 
tron spin magnetization orientated along the direction 
of the static magnetic field (1/r = 1/T\ e + 1/T2 e ), T\ e 
is the electron spin lattice relaxation time and T2 e the 
transverse electron spin relaxation time. The distance 
f3 = (C/D) 1 / 4 (D is the spin-diffusion constant), charac- 
terizes the cross-over radius between direct and diffusion- 
mediated hyperpolarization [37] . 

For nuclei located outside this radius (r > j3) , polariza- 
tion occurs predominantly via nuclear spin diffusion with 
a rate r^ 1 = ^ = DV 2 P [14J. D = Wo 2 ~ a 2 /(50T 2n ) 
is the spin-diffusion constant, with a the average sepa- 
ration between nearest-neighbor nuclei, W the probabil- 
ity of a flip-flop transition between nuclei due to dipole- 
dipole interaction, and T2 n the transverse nuclear spin 
relaxation time [37J. For the present system, we esti- 
mate — 3 nm and /?si — 4 nm, given T\ e c± 30/is 
and T 2e ~ 1 ns [3J, a H (si) ~ 1-3 (4.1) A, and T 2nH (Si) = 
0.01 (5.6) ms. Transmission electron micrograph studies 
of the SiNP show an oxide thickness of a few nm, a dis- 
tance small enough for this direct polarization process to 
take place. We note that j3 is not a hard cutoff of the 
distance for direct polarization to take place, rather it 
gives a transition point between polarization dominated 
by direct DNP and spin diffusion. 

Although the total enhancement (ranging from ~ 1.5 — 
3) over all X H spins in the sample is relatively low com- 
pared to enhancements in concentrated radical solutions 
[38J, the enhancement near the surface of the particles 
is significant. Numerical modeling of the polarization as 
a function of distance from the particle surface shows a 
characteristic length of order 10 nm. This is consistent 
with studies of polarization as a function of radical con- 
centration [T5]. Estimating the enhanced polarization of 
X H nuclei within 10 nm of the surface of the SiNP (cor- 
responding to ~ 1.2% of all X H nuclei in the sample) 
based on the model above, we find a factor of ~ 50 above 
equilibrium polarization at 4.2 K for the natural abun- 
dance particles in the ethanol solution, corresponding to 
an enhancement of ~ 3000 times the room temperature 
equilibrium polarization. This enhancement exceeds the 
reported enhancements of functional groups on silica sur- 
faces using artificial free radicals and much higher power 
microwave sources, estimated at ~ 30 times the equilib- 
rium polarization at 77 K for molecules within 1 nm of 
the surface [39J. 

This effect may be optimized for enhanced solution po- 
larization or studies of surface functionalization by reduc- 
ing the 29 Si concentration, increasing the surface area of 
the SiNP, using porous silicon, or engineering thinner 
surface oxides with higher defect densities. Additionally 
we note that this effect may be revealed in other solid 
state systems with a lower nuclear spin concentration or 
lower gyromagnetic ratio. In particular, 13 C DNP has 



been demonstrated in a variety of diamond materials [2 J 
(1% 13 C natural abundance) and nanodiamond and has 
surface active electronic defects |40] that may be suit- 
able as polarizing agents by mechanisms similar to those 
described in this paper. 
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